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ABSTRACT 


The  purpose  of  this  thesis  is  to  design  a  quantum  well  infrared  detector  with 
tunable  spectral  bandwidth.  The  tunability  of  the  bandwidth  is  achieved  by  using  the 
linear  Stark  effect  for  the  ground  to  first  excited  state  transition  in  an  asymmetric 
quantum  well.  The  position  of  the  absorption  peak  is  dependent  on  the  direction  of  the 
electric  field,  and  therefore  it  can  be  either  blue  or  red  shifted  by  changing  the  direction 
of  the  field.  If  two  identical  asymmetric  quantum  wells  are  arranged  opposite  each  other, 
we  can  obtain  both  the  blue  and  red  shift  for  either  direction  of  the  bias.  This  method  can 
produce  broader  peaks  with  tunable  bandwidths  proportional  to  the  applied  field. 

A  program  was  developed  to  calculate  the  energy  levels  and  wavefunctions  of  an 
arbitrary  quantum  well.  The  program  was  used  to  design  a  step  quantum  well  capable  of 
detecting  infrared  in  the  8-12  [im  band.  The  validity  of  the  approach  was  verified  by 
comparison  with  experimental  data  of  Martinet  et  al.  (Martinet,  1992),  and  found  to  have 
a  good  agreement.  The  designed  step  well  was  used  to  create  a  tunable  bandwidth 
detector.  The  analysis  showed  that  the  bandwidth  could  be  tuned  to  more  than  twice  the 
peak  width.  The  numerical  simulation  indicates  the  possibility  of  manufacturing  a  tunable 
bandwidth  infrared  detector  by  using  step  quantum  wells. 
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I.  INTRODUCTION 


A.  INFRARED  DETECTORS 

The  human  eye  can  detect  only  a  small  segment  of  the  electromagnetic  spectrum, 
the  familiar  visible  band  ranging  from  blue  to  red.  Only  objects  hot  enough  to  radiate  at 
these  frequencies  and  those  that  reflect  visible  light  are  visible  to  the  eye.  For  the  most 
part,  however,  objects  are  too  cold  to  glow  noticeably,  making  them  undetectable  at 
night.  Yet  their  finite  temperature  does  give  them  an  infrared  glow.  It  is  this  glow  that 
night- vision  infrared  cameras  strive  to  see.  At  the  heart  of  an  infrared  camera  is  its  “eye”: 
a  two-dimensional  array  of  detector  pixels,  with  each  pixel  converting  some  of  the 
incident  photons  to  an  electric  signal. 

A  Quantum  Well  Infrared  Photodetector  (QWIP)  consists  of  a  stack  of  quantum 
wells  designed  to  detect  infrared  light  (Levine,  1993).  When  the  quantum  well  is 
adequately  deep  and  narrow,  its  energy  states  are  quantized.  The  potential  depth  and 
width  of  the  well  can  be  adjusted  so  that  it  has  only  two  energy  states:  a  ground  state  near 
the  well  bottom  and  a  first  excited  state  near  the  well  top.  In  Figure  1.1,  photoexcited 
electron  transition  from  the  ground  state  to  the  first  excited  state,  where  an  externally 
applied  voltage  sweeps  them  out,  producing  a  photocurrent. 


Figure  1.1  Electron  excitation  of  a  quantum  well  infrared  photodetector. 

For  two  bound  states  in  the  quantum  well,  only  photons  having  energies 
corresponding  to  the  energy  separation  between  the  two  states  are  absorbed,  resulting  in  a 
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detector  with  a  sharp  absorption  spectrum.  If  there  is  only  one  bound  state  in  the  well,  the 
transition  occurs  between  the  bound  state  and  the  continuum,  producing  a  broader  peak 
(Levine,  1993).  Designing  a  quantum  well  to  detect  light  of  a  particular  wavelength 
becomes  a  straightforward  matter  of  adjusting  the  potential  depth  and  width  of  the  well  to 
produce  two  states,  separated  by  the  desired  photon  energy.  Impressive  advances  in  the 
crystal  growth  and  processing  techniques  of  compound  semiconductors  such  as  gallium 
arsenide  (GaAs)  and  aluminum  gallium  arsenide  {Alfia^_^As)  with  different  aluminum 

compositions  allowed  the  fabrication  of  novel  QWIP  detectors  (Levine,  1993)  and  arrays 
(Gunapala,  1996). 

B.  PROGRESS  OF  QWIPS 

In  the  1970's,  InSb  and  HgCdTe  were  used  to  cover  the  3-5  and  8-12  micrometer 
wavelength  ranges.  However,  HgCdTe  has  some  shortcomings  due  to  difficulty  in  its 
material  preparations,  uniformity  and  reproducibility  in  the  fabrication  process.  Other 
difficulties  include  mechanical  softness,  sensitivity  to  elevated  temperature,  and  low  yield 
for  large  array  fabrications,  which  make  such  detectors  extremely  expensive.  In  order  to 
fabricate  low-cost  detectors  and  overcome  the  above  difficulties,  detection  using 
superlattices  was  suggested  by  Esaki  and  Tsu  (Esaki,  1969).  Extensive  investigations, 
both  theoretical  (Coon,  1984  and  1985)  and  experimental  (Smith,  1983),  were  a  subject 
of  considerable  interest  and  were  carried  out  until  the  first  observation  of  strong 
intersubband  transitions  in  GaAsI  AlGaAs  quantum  wells  by  West  and  Englash  (West, 
1985).  Since  the  problems  associated  with  material  preparation  were  avoided,  this  new 
development  has  received  considerable  attention. 

The  first  bound-to-bound  quantum  well  infrared  photodetector  was  demonstrated 
by  Eevine  et  al.  (Eevine,  1987).  Because  the  photoelectrons  had  to  tunnel  out  of  wells 
from  the  bound  excited  state  under  large  barriers,  as  illustrated  in  Eigure  L2(a),  the 
responsivity  was  so  low  that  the  device  had  to  work  using  a  CO2  laser.  This  low 
responsivity  led  to  the  use  of  bound -to- continuum  transitions,  which  were  first  proposed 
by  Coon  and  Karunasiri  in  1984  (Coon,  1984).  Hasnain  et  al.  (Hasnain,  1989) 
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demonstrated  the  first  bound -to- continuum  QWIP  (see  Figure  1.2b)  to  have  a  much 
higher  responsivity,  but  also  a  small  peak  absorbance  due  to  broad  spectral  width. 


Figure  1.2  Schematic  diagrams  showing  (a)  bound -to- bound  and  (b)  bound -to- 

continuum  electron  transfer. 

Several  configurations  of  the  QWIP  structure  have  been  reported  based  on 
transitions  from  bound -to- extended  states  (Levine,  1990),  bound-to-quasi-continuum 
states  (Levine,  1991),  bound-to-quasibound  states  (Gunapala,  1996),  and  bound-to- 
miniband  states  (Yu,  1991).  These  configurations  provide  good  detection  performance 
because  in  all  designs,  electrons  in  the  ground  state  do  not  flow  in  response  to  the  bias, 
while  electrons  activated  to  the  excited  state  yield  photocurrent  under  low  bias.  Several 
structural  optimizations  have  also  been  made  to  improve  the  performance  of  QWIPs  such 
as  dark  current  reduction  from  thermionic  emission  (Gunapala,  1996). 

C.  PURPOSE  OF  THIS  THESIS 

The  goal  of  this  research  is  to  design  a  quantum  well  detector  with  a  tunable 
spectral  bandwidth.  The  tunability  of  the  bandwidth  is  achieved  by  using  “the  linear  Stark 
effect”  in  an  asymmetric  quantum  well  (Yuh,  1989).  For  example,  in  the  step  quantum 
well  shown  in  Figure  L3(a),  the  amount  of  Stark  shift  from  the  ground  to  first  excited 

bFu 

states  transition  is  approximately  ,  where  F  is  the  electric  field  strength,  and  a  is 
the  step  width  (Yuh  1989). 
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Figure  1.3  Schematic  diagrams  showing  an  (a)  asymmetric  step  quantum  well  and  a 
(b)  designed  step  quantum  well  structure  for  bandwidth  tuning. 

As  previously  mentioned,  the  position  of  the  absorption  peak  is  sensitive  to  the 
direction  of  the  electric  field  and  can  be  either  blue  or  red  shifted  by  changing  the 
direction  of  the  field  (Mii,  1988).  If  two  wells  are  arranged  as  in  Figure  1.3(b),  both  red 
and  blue  shifts  can  occur  at  the  same  time  producing  an  absorption  peak  with  a  tunable 
bandwidth. 

In  the  following  chapter,  we  will  present  a  detailed  analysis  of  energy  levels  in  a 
step  quantum  well,  including  external  electric  field  effects.  In  addition,  the  absorption 
coefficient  under  applied  field  will  be  calculated  using  a  step  approximation  to  the  linear 
potential.  Finally,  in  Chapter  3,  we  will  present  an  optimized  device  structure  to  achieve 
tunable  bandwidth  detection. 
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11.  THEORETICAL  BACKGROUND 


A.  INTRODUCTION 

In  the  following  sections,  we  will  develop  the  theoretical  tools  necessary  to 
analyze  quantum  well  energy  levels  as  a  function  of  applied  bias  for  an  asymmetric 
quantum  well.  The  Schrodinger  equation  is  solved  numerically  using  the  transfer  matrix 
approach  (Wang,  1989).  The  effect  of  the  electric  field  on  the  energy  levels  will  be 
incorporated  by  the  “staircase”  approximation  of  the  linear  potential. 


B.  TRANSFER  MATRIX  METHOD 

In  the  case  of  an  infinite  square  quantum  well,  the  Schrodinger  equation  is  given 

by 


= 


n  =  1,2,3,... 


(2.1) 


Im^dz 

where  zis  the  growth  direction,  n  is  the  quantum  number,  ti  is  the  reduced  Planck’s 
constant,  is  the  effective  mass  of  the  electron,  and  and  are  the  wavefunction 


and  energy  in  the  energy  level,  respectively.  The  wave  function  ?  must  also  satisfy 

boundary  conditions  across  a  potential  step  since  the  probability  current  across  a 
boundary  should  be  continuous. 


In  light  of  the  above  the  conditions  and  the  fact  that  the  wave  functions  go  to  zero 
at  the  boundaries  (z  =  0,L)  and  the  electron  stays  within  the  well,  we  can  obtain  the 
energy  and  wavefunction  of  the  quantized  energy  states  as 

n  =  1,2,3,...  (2.2) 

2m,  L 

and 


^^1^)  =  Asin 


n  =  1,2,3,... 


(2.3) 
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where  L  is  the  width  of  the  quantum  well.  Taking  the  origin  at  the  left  side  of  the  well  as 
shown  in  Figure  2.1,  the  wave  function  can  then  be  normalized  over  the  total  distance 
-  oo  <  z  <  +00  as  follows: 

j'F:(z)'F„(z)dz  =  l,  (2.4) 


Figure  2.1  Infinite  square  well  potential  and  energy  levels. 

There  are  many  techniques  available  to  calculate  the  energy  levels  and  the  wave- 
functions  of  finite  quantum  well  structures,  each  having  advantages  and  disadvantages. 
The  deciding  factor  on  which  method  to  use  depends  on,  to  some  extent,  the  structure 
itself.  For  example,  the  variational  approach  yields  simple  results  for  the  ground  state  but 
has  difficulty  finding  the  excited  states.  The  solution  of  the  Schrodinger  equation  for  a 
finite  square  quantum  well  can  be  found  by  matching  boundary  conditions  (Weisbuch, 
1991).  However,  when  the  number  of  interfaces  is  large,  this  method  becomes 
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cumbersome  (Choi,  1997).  With  the  help  of  a  computer,  the  transfer- matrix  method  is  a 
multipurpose  approach  in  obtaining  the  energy  states  and  the  corresponding  wave 
functions  of  a  structure  with  an  arbitrary  potential  profile.  In  the  following  section,  the 
transfer  matrix  method  will  be  discussed. 

A  heterostructure  is  made  from  layers  of  materials  having  different  bandgaps.  For 
example,  in  a  typical  square  quantum  well,  a  layer  of  GaAs  is  surrounded  between  the 
two  thicker  higher  bandgap  layers  of  Alfia^_^As .  The  well  width  is  controlled  by 
adjusting  the  GaAs  layer  thickness  while  the  potential  depth  is  controlled  by  adjusting 
the  Al  composition  in  the  barrier  layers.  Again  the  application  of  the  boundary 
conditions  leads  to  the  quantized  energy  values  and  wave  functions  . 


In  the  transfer  matrix  method,  a  plane  wave  is  assumed  to  be  incident  from  the  left 
of  the  structure  as  shown  in  Figure  2.2.  In  each  layer  the  wave  function  can  be  written  as 

a  linear  combination  of  the  forward  ( )  and  backward  ( )  traveling  plane 

waves.  A  matrix  can  be  set  up  using  the  boundary  conditions  at  each  interface  of  the 
quantum  weU  structure  to  connect  the  wavefunctions  between  the  adjacent  regions 
(Hutchings,  1989). 
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Figure  2.2  Arbitrary  multi-layered  potential  well  structure. 

For  the  arbitrary  quantum  well  structure  shown  in  Figure  2.2,  the  wave  function 
for  the  n'*  region  can  be  written  as 
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n  =  \,2, . ,N 


(2.6) 


where  A„and  are  arbitrary  constants, 


K=J^(E-V„).  (2.7) 

m*  is  the  electron  effective  mass  in  the  n'*  layer  and  the  thickness  of  the  layer 
-z„.  Furthermore,  E  is  the  electron  energy,  which  can  either  be  higher  or 

lower  than  the  highest  potential  in  the  structure.  The  boundary  conditions  at  the  n‘*' 
interface  give  (Robinett,  1997) 

¥„(4)=¥„+i(zJ 


1  d\]f„{z„)_  1  d\^„^i(4) 

K  dz„  dz„ 

Thus,  the  following  matrix,  ,  can  be  obtained  to  connect  the  coefficients  in  the  n* 
layer  to  that  of  the  (n  + 1)'*  layer  as  follows: 


where 


and 


,1 

1 

n+l 

=  M„ 

n 

_^n+l  _ 

A. 

n  =  1,2,... ,n  -1 


M„  =- 
"  2 


(l  +  y  likn+kj^+])Zn 


I  n  J  * 


(2.9) 


(2.10) 


(2.11) 


In  Equation  (2.10)  is  the  spatial  coordinate  of  the  interface  between  n'*  and  (n  +  1)'' 


layers. 


This  model  can  be  extended  to  the  multi-layered  quantum  well  structure  and  by 
using  the  above  relationship,  the  coefficients  of  the  wavefunctions  for  the  first  and  last 
layers  are  related  by 


A 

=  M 

'A' 

_ 

mil 

'A' 

A. 

_m2i 

"^2. 

A. 

(2.12) 
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If  the  two  outer  layers  have  higher  potentials  than  the  energy  of  the  particle  (i.e.,  if 
E  <  V  ),  then  there  will  be  bound  energy  states.  In  this  case,  we  can  set  the  coefficient 
Bff  =  0 ,  because  it  would  become  infinite  as  z  — >  oo .  In  the  same  manner,  A^  can  be  set 
to  zero  for  the  region  z  <  0  .  Under  these  conditions.  Equation  2.12  reduces  to 

(2.13) 


(2.14) 


A.' 

m,i 

m,2 

'o' 

0 

_m2i 

A. 

In  order  to  satisfy  the  Equation  2.13,  we  must  have 

m^^iE)  =  0. 


A  Matlab  program  was  developed  to  plot  (E)  versus  E  and  identify  the  points  where 
m22  is  zero.  Eor  example,  Eigure  2.3  shows  the  calculated  111^2  as  a  function  of  E  for  a 
square  potential  well.  The  calculated  bound  state  energies  and  wavefunctions  are  shown 
in  Eigure  2.4. 


-IQ*'  Finding  Energy  Levels  In  a  Pocential  Wei 


0 

Eigure  2.3  Plot  of  m22(E)  versus  E  for  a  square  potential  well  with  L  =  110  A, 

Vj  =  =  0.75  eU  and  y2  =  0  eU  . 
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Figure  2.4  Energy  levels  and  wavefunctions  of  the  quantum  well  described 

in  Figure  2.3. 


C.  INTERSUBBAND  TRANSITIONS 

In  a  quantum  well  structure  made  of  semiconductors,  the  wavefunction  consists  of 
both  envelope  and  Bloch  functions.  Thus,  not  only  the  transitions  between  Bloch 
functions  but  also  between  the  envelope  states  are  allowed  within  the  structure.  The 
transitions  between  Bloch  states  correspond  to  interband  transitions,  while  transitions 
between  envelope  states  correspond  to  intersubband  transitions.  For  a  given  interaction 
potential  the  transition  rate  W  can  be  obtained  from  the  Fermi’s  golden  rule  (Coop 

1984) 

W  =  (2.15) 

^  / 


10 


where  and  are  the  electron  energies  in  the  initial  and  final  states,  respectively,  and 
/zoo  is  the  incident  photon  energy.  The  interaction  potential  energy,  V^,  in  the  dipole 
approximation  is  given  by  (Coon,  1984) 


in 


xl/2 


2eQ«^coc 


e-  p 


(2.16) 


where  e  is  the  electron  charge,  /  is  the  incident  photon  flux,  c  is  the  speed  of  light,  8o 
is  the  electric  permittivity,  is  the  index  of  refraction,  CO  is  the  angular  frequency  of  the 


incident  photon,  e  is  the  unit  vector  of  the  photon  polarization  direction,  and  p  is  the 
momentum  operator. 


p  =  -in 


6^  d  ^  d  ^  d  ^ 

i - 1-  / - \-k  — 

^  dx  dy  dz  j 


(2.17) 


In  the  case  of  a  quantum  well,  the  wavefunctions  for  the  initial  and  final  states  of 
transition  can  be  written  as 


(2.18) 


where  M,and  M,are  the  Bloch  functions  of  the  initial  and  final  states,  respectively. 


k^  and  /r.  are  the  initial  and  final  wave  vectors  in  the  xy  plane,  r  is  the  3-D  position 


vector,  p  is  the  position  vector  in  the  xy  plane,  'B,  and  are  the  corresponding 

envelope  functions,  and  A  is  the  area  of  the  structure.  The  transition  matrix  element  in 
Equation  2.15  can  be  evaluated  as  (Coon,  1984) 
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< I I 'P,  >=<  I  I  M,(r>|/,(z)^ 


^<u  (r)e''^''  ly  IM-Cry’^'p  ><\|/  (z)l\|/.(z)>  (2.19) 


+  <\|/^(z)^^iypl\|/.(z)^><  M^(r)lM.(r)>. 


In  the  above  expression,  we  assume  that  the  envelope  funetions  vary  slowly  eompared  to 
the  Bloeh  functions.  The  integration  over  the  x  and  y  coordinates  gives  k-  =  kj^  .  The  first 
term  of  Equation  2.19  is  zero  because  of  the  fact  that  the  Hamiltonian  used  to  obtain  \|/, 
and  V|/y  is  Hermitian.  In  addition,  the  Bloch  functions  are  normalized  to  unity.  Thus,  the 
intersubband  transition  matrix  element  reduces  to  (Coon,  1984) 


<'y^iy^i'y,>=<vi//z)iy,iv|/,(z)> 


in 


-sl/2 


(2.20) 


<\|/.(z)IPJ\|/,(z)>cose, 


where  0  is  the  angle  between  the  photon  polarization  angle  and  the  z  axis  and  P  is  the 


momentum  operator  in  the  z  direction.  The  strength  of  the  transition  is  usually 
characterized  by  the  oscillator  strength,  which  is  defined  as 


where 


/  = 


r 


(2.21) 


('E^lzl'E,)  =j  'P^z'E,dz. 


(2.22) 


D.  INTERSUBBAND  TRANSITION  SELECTION  RULES 

The  first  selection  rule  of  intersubband  transition  is  the  parity  selection  rule.  From 
Equation  2.20,  it  is  found  that  the  parity  of  the  final  state  has  to  be  different  from  that  of 
the  initial  state  in  order  to  obtain  a  non-  zero  transition  matrix  element.  This  means  that 


12 


the  intersubband  transitions  in  a  square  quantum  well  can  occur  only  from  odd  states  to 
even  states  or  vice  versa.  For  example  1  ^  2  or  2  — >  3  transitions  are  allowed  whereas 
1— >3  is  not  allowed.  Figure  2.5  shows  schematically  the  parity  selection  rule  of  the 
intersubband  transitions  in  a  square  potential  well.  However,  the  symmetry  of  a  square 
well  can  be  broken  by  applying  a  relatively  strong  electric  field,  making  the  normally 
“forbidden”  transitions  become  allowed  (Mii,  1990). 
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Figure  2.5  Square  quantum  well  showing  allowed  and  forbidden  transitions. 


The  second  selection  rule  is  the  polarization  selection  rule.  Since  the  potential 
modulation  of  a  quantum  well  structure  is  only  in  one  direction,  the  z  direction, 
intersubband  transition  will  not  occur  when  the  photon  polarization  is  perpendicular  to 
the  z  direction.  This  can  be  understood  by  considering  the  cos0  term  in  Equation  2.20. 
This  unique  polarization  selection  rule  is  an  important  characteristic  and  can  be  used  for 
the  identification  of  the  intersubband  transition  (West,  1985).  In  order  to  obtain  a  strong 
intersubband  transition,  the  photon  polarization  has  to  coincide  with  the  growth  direction 
(i.e.,  0=0). 

For  an  asymmetric  quantum  well  (see  Figure  2.6),  the  dipole  matrix  elements  of 
transitions  for  both  odd-to-odd  and  even-to-even  quantum  numbers  do  not  vanish,  since 
the  eigenfunctions  of  these  energy  states  are  not  eigenfunctions  of  the  parity  operator. 
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Figure  2.6  Asymmetric  quantum  well  with  allowed  transitions. 


In  addition  to  allowing  intersubband  transitions  between  any  low  states,  the 
asymmetric  well  structure  also  gives  linear  Stark  shift  (Yuh  1989).  This  can  be  utilized 
for  tuning  the  detection  wavelength  (Martinet ,  1992). 


E.  INTERSUBBAND  ABSORPTION  COEFFICIENT 

The  absorption  coefficient  in  quantum  well  devices  is  one  of  the  most  important 
factor  in  the  design  of  the  quantum  structures  and  is  proportional  to  the  amount  of 
incident  light  that  is  absorbed  in  the  structure.  The  absorption  coefficient  depends  upon 
the  material  characteristics  of  the  respective  semiconductor.  For  example,  the  high 
refractive  index  of  a  semiconductor  such  as  GaAs ,  (n  =  3.3),  causes  that  light  incident  on 
the  surface  of  an  epitaxial  layer  structure  is  refracted  towards  the  normal,  for  example, 
and  away  from  the  required  orientation.  Efficient  coupling  of  the  intersubband  transition 
can  be  achieved  by  a  number  of  optical  techniques,  including  the  use  of  a  beveled  edged 
device  (Levine,  1987),  diffraction  gratings  (Hasnain,  1989),  and  combined  diffraction 
gratings  and  waveguides  (Andersson,  1992). 

In  addition  to  absorption,  it  is  important  to  extract  the  photoexcited  carriers  out  of 
the  quantum  well.  This  is  achieved  by  placing  the  excited  state  near  the  top  of  the  barrier 
as  shown  in  Figure  2.7.  In  the  case  of  the  bound-to-continuum  transition  (see  Figure 
2.7a),  the  excited  state  is  located  just  above  the  barriers,  while  for  a  bound-to-quasi 
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bound  structure  (see  Figure  2.7b)  the  excited  state  is  just  below  the  barrier.  The  latter 
structure  provides  lower  leakage  current  due  to  the  higher  barrier  height  for  the  same 
detection  wavelength. 


(a) 


(b) 


Figure  2.7  Bound-to-continuum  (a)  and  (b)  bound -to- bound  transitions. 

For  the  bound -to- bound  transitions,  since  there  is  only  one  final  state,  the 
summation  over  the  initial  states  in  Equation  2.23  gives  the  number  of  the  electrons  in  the 
ground  state. 


N  =  N^LA .  (2.23) 

Here,  N  is  the  number  of  electrons  in  the  initial  state,  is  the  three  dimensional 

doping  density,  and  L  and  A  are  the  well  width  and  area  of  the  quantum  structure, 
respectively.  The  delta  function  in  Equation  2.15  is  usually  replaced  by  the  normalized 
Eorentzian  function  to  take  into  account  the  spectral  line  width  given  by  (Yariv,  1997) 


g(h(0)  = 


J _ T _ 

2"  (£,-£, 


(2.24) 
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where  T  is  the  full  width  at  half  the  maximum.  This  is  usually  obtained  experimentally 
and  is  about  10-20  meV  for  intersubband  absorption.  The  transition  rate,  found  using  the 
interaction  potential  and  the  momentum  operator  is  given  by 


'  Bound  Bound 


N,{IAL)e^h^ 


k^J  — IV,>Pg(;i(0). 


(2.25) 


For  the  bound -to- continuum  transitions,  the  number  of  initial  states  is  still  N 
with  energy  ,  as  in  the  case  of  bound -to- bound  states.  However,  the  summation  over 
the  final  states  requires  multiplication  by  the  joint  density  of  states  p  {E) ,  which  is 


defined  as 


p{E)dE- 


(2.26) 


For  energy  states  separated  by  the  incident  photon  energy,  the  above  equation  can  be 


reduced  to 


2iiih\E,-V,) 


(2.27) 


In  this  expression,  I  is  the  length  of  the  region  used  for  the  normalization  of  the  extended 
state  wavefunction,  ml  is  the  electron  effective  mass  in  the  barrier,  and  I4  the 
potential  height  of  the  barrier.  The  transition  rate  for  bound -to- continuum  states  is  then 


defined  as 


'  Bound  Continuum 


N/IAL)e\^l  2m*  d  2 

^ - - -  \<WA - l\|/.  >r . 

2(m:)  dz 


(2.28) 


This  model  of  the  bound  to  free  intersubband  absorption  is  one  where  the  line  width  is 
due  only  to  the  variation  of  the  joint  density  of  states.  The  absorption  coefficient  a(^(0) 
can  be  obtained  by  using  both  the  transition  of  the  matrix  element  and  the  transition  rate 
as  follows: 


a(fiQd)  = 


number  of  transitions  per  unit  volume  and  time  ,  /at 

- ^ ^ - - - ,  a(h(a)= 


(2.29) 


incident  photon  flux  / 

Thus,  the  absorption  coefficients  for  the  bound -to- bound  and  bound -to- continuum 
transitions  are  given  by 
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2fc3 


a 


Bound— ¥  Bound 


(hCd)  = 


|<V|/  \-^\\\f.  >P  gifKO) 

dz 


(2.30) 


a 


Bound  —^Continuum 


(fidi)  = 


N,eWl 


2m^ 


2(m^)  8on^c 


l<\|/ ,  I  — 1\|/,  >P 


The  length  term,  /,  which  appears  in  the  bound -to- continuum  transition,  is  eliminated  by 
the  ^  term  in  the  normalized  final  state  wave  function  (assuming  it  is  confined  to  a 


length,  /).  As  an  example,  absorption  coefficients  for  the  bound-to-bound  and  bound-to- 
continuum  situations  for  the  quantum  well  shown  in  Figure  2.8  are  depicted  in  Figures 
2.9  and  2.10,  respectively. 
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Figure  2.8  Energy  levels  of  the  AIq^Ouq^^As  /  / GaAs  /  Al^^Ga^f^As 

quantum  potential  well. 
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Figure  2.9  Bound -to- bound  absorption  coefficient  for  the  quantum  potential  well 

given  in  Figure  2.8  with  Y  =  20  meV  . 


Figure  2.9  shows  the  absorption  versus  photon  energy  for  bound -to- bound 
transitions.  There  is  a  sharp  rise  in  the  absorption  spectrum  when  the  photon  energy  is 
equal  to  the  difference  between  the  bound  state  energies  in  the  quantum  well  structure. 
The  width  of  the  curve  depends  on  the  finite  lifetime  of  the  excited  electron,  which  was 
modeled  using  the  Gaussian  spectral  density  function. 
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Bowid  to  Continuum  Absorption 


Figure  2.10  Bound- to- continuum  absorption  coefficient  for  the  quantum  well 

in  Figure  2.8. 

Similarly,  Figure  2.10  shows  the  absorption  coefficient  versus  the  photon  energy 
for  the  case  of  bound -to- continuum  transitions.  The  sharp  rise  in  the  spectrum 
corresponds  to  a  photon  energy  equal  to  the  difference  between  the  bound  state  energy 
and  the  energy  at  the  top  of  the  quantum  well  where  the  extended  continuum  states  begin. 
As  the  photon  energy  increases  above  this  threshold  value,  absorption  increases  until 
reaching  a  peak  value  corresponding  to  the  energy  above  the  top  of  the  barrier.  At  very 
high  energies  there  is  a  slow  decrease  of  a  (/zoo)  since  the  density  of  states  drops  as  the 
energy  is  increased. 

As  previously  mentioned,  our  goal  is  to  design  a  detector  with  tunable  bandwidth. 
The  most  convenient  way  to  do  this  is  by  using  bias  as  the  tuning  parameter.  As  a  result, 
it  is  important  to  study  the  effect  of  bias  on  the  energy  levels  as  well  as  on  the  absorption. 
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F.  EFFECT  OF  ELECTRIC  FIELD  ON  QUANTUM  WELLS 

When  a  uniform  electric  field  is  applied  along  the  z  direction  (growth  direction), 
the  quantum  well  structure  will  tilt  in  the  opposite  direction  of  the  electric  field  and  the 
change  due  to  this  tilt  will  affect  the  energy  states  of  the  quantum  well.  The  potential 
experienced  by  an  electron  in  the  quantum  well  is  given  by 

V{z)=V„  +  eFz,  (2.31) 

where  is  the  potential  Mght  of  the  n*  layer  of  the  quantum  well,  e  is  the  electron’s 
charge,  F  ,  measured  in  (U  /m) ,  is  the  electric  field  strength  ,  which  is  applied  along  the 
growth  direction.  A  square  quantum  well  in  the  presence  of  an  external  bias  is  shown  in 
Figure  2.11. 


Figure  2. 1 1  Schematic  representation  of  a  square  quantum  well  under  an  applied  bias. 

By  increasing  the  electric  field,  the  lowest  energy  state  in  the  quantum  well 
will  decrease  in  energy  and  fall  into  to  the  triangular  potential  formed  at  the  bottom  of  the 
well.  On  the  other  hand,  the  excited  states  are  less  affected  by  the  external  electric  field. 
However,  for  low  electric  field  strengths,  the  shifting  of  all  the  energy  levels  are  the  same 
due  to  the  symmetry  of  the  square  well.  Thus,  there  is  no  appreciable  Stark  effect.  This 
can  be  altered  using  an  asymmetric  quantum  well,  which  is  equivalent  to  a  square  well 
under  a  high  electric  field. 
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1.  Transfer  Matrix  Approach  for  Quantum  Wells  Under  Applied  Bias 

Previously,  we  have  discussed  the  calculation  of  the  energy  levels  and 
wavefunctions  of  arbitrary  quantum  wells  using  transfer  matrix  approach.  However,  this 
approach  cannot  be  used  directly  due  to  the  linear  variation  of  the  potential  energy  under 
an  external  bias. 

One  way  to  overcome  this  problem  is  to  approximate  the  linear  potential  due  to 
the  bias  shown  in  Figure  2.12,  with  a  series  of  steps  or  “staircase”  as  shown  in  Figure 
2.13.  The  advantage  of  this  approach  compared  to  the  use  of  Airy  functions  is  that,  it  can 
be  used  at  low  field  strengths,  whereas  Airy  functions  tend  to  diverge  at  low  bias 
voltages,  making  the  numerical  analysis  much  more  difficult  (Brennan,  1987). 


Quantum  Wd  with  5x10®  VAn  Extefnal  Electric  FieW 


Figure  2. 12  Step  quantum  well  with  applied  bias  of  5  F  /  pm . 
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Figure  2.13  The  stair  case  approximation  of  the  potential  shown  in  Figure  2.12. 

The  Schrodinger  equation  for  an  electron  in  the  quantum  well  in  the  presence  of 
an  electric  field  F  can  be  written  as 

-^^'V,-{E-V,-qFz)'V,=0,  (2.32) 

2m„  dz 

where  E  is  the  energy  of  the  electron  in  the  layer  and  and  are  the  potential  and  the 


effective  mass  in  the  n*  layer,  respectively.  In  order  to  apply  the  transfer  matrix  relations 
it  is  necessary  to  know  the  boundary  conditions  at  the  end  points.  For  bound  states,  we 
can  use  the  same  approach  as  used  in  the  zero  bias  case.  For  the  continuum  states  we  can 
assume  that  the  electron  is  incident  from  the  left  side  of  the  potential  shown  in  Figure 


2.13,  and  which  gives  5^=0.  Thus,  starting  from 


all  the  unknown 
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coefficients  can  be  progressively  calculated,  ending  with  the  coefficient 

N 

.4_ 

0 

Therefore  the  electron  wave  functions  can  easily  be  obtained  because  all  the  coefficients 
are  produced  when  is  known  (Hutchings,  1989). 
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Figures  2.14  (a)  through  2.14  (d)  show  the  potentials  and  eigenenergies  ealeulated 
using  the  “stairease”  approximation  for  external  bias  varied  from  0  to  5F/|lm.  In  this 
way,  we  ean  observe  the  ehange  in  the  energy  states  as  a  result  of  the  applied  eleetrie 
field.  Table  2.1  shows  the  ehange  in  bound  state  energies  due  to  the  external  applied  bias 
in  the  quantum  well. 
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Figure  2.14  Effeet  of  electric  field  on  the  bound  state  energies  of  a  step  quantum  well 
for  the  following  field  strengths:  (a)  F  =  0 ,  (b)  F  =  1  V /  \im ,  (c)  F  =  3  VI  pm  and 

(d)  F=  5  F/pm. 
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Zero 

IE/  |lm 

3E/  |lm 

5V  !  |lm 

Energy  States 

Biased  Well 

Biased  Well 

Biased  Well 

Biased  Well 

(meV) 

(meV) 

(meV) 

(meV) 

47 

61 

74 

El 

41 

?E=6 

?E=20 

?E=33 

153 

164 

175 

E2 

147 

?E=6 

?E=17 

?E=28 

250 

256 

-- 

E3 

247 

?E=3 

?E=9 

Continuum 

Table  2.1  Bound  state  energies  as  a  function  of  bias. 


2.  Intersubband  Absorption  Under  an  Applied  Electric  Field 


The  intersubband  absorption  peak  position  will  shift  in  energy  with  changing 
electric  field  due  to  the  shift  of  the  energy  levels.  At  a  low  bias,  this  effect  is  especially 
pronounced  in  asymmetric  quantum  wells  as  compared  to  symmetrical  wells.  Figure  2.15 
shows  the  calculated  absorption  coefficient  as  a  function  of  bias  across  the  well  using 
Equation  2.30.  The  matrix  element  in  Equation  2.30  is  evaluated  using  the  wavefunction 
under  the  applied  bias. 
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3000 


Boutx)  to  Bouid  Absorpton  With  Applied  Bias 


Figure  2.15  Absorption  coefficient  as  a  function  of  bias  across  the  step  quantum  well 

structure. 

The  strength  of  absorption  is  found  to  increase  with  increasing  electric  field, 
mainly  due  to  the  stronger  matrix  element.  The  asymmetry  observed  between  positive 
and  negative  biases  is  due  to  the  fact  that  under  forward  bias  (negative),  the  first  excited 
energy  level  moves  toward  the  ground  state  and  hence  the  transition  energy  is  reduced, 
while  under  reverse  bias  (positive)  the  opposite  is  true  and  the  transition  energy  is 
increased. 

G.  COMPARISON  WITH  EXPERIMENT 

In  order  to  verify  the  validity  of  the  “staircase”  model  for  a  quantum  well  under 
bias,  we  have  analyzed  a  step  quantum  well  structure  employed  by  Martinet  et  al. 
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(Martinet  1992)  for  the  study  of  bias  dependent  photocurrent.  The  parameters  of  the 
structure  used  in  Martinet’s  experiment  are  shown  in  the  Table  2.2. 


Entered  Experimental  Data 

Layer 

and  Thickness 

Concentration 

AlGaAs  (barrier) 
36.5  nm 

0.40 

AlGaAs  (step)  6.7 

0.20 

GaAs 

3.2nm 

— 

AlGaAs  (barrier) 
36.5  nm 

0.40 

Table  2.2  Parameters  of  the  step  quantum  well  structure  used  by  Martinet  et  al. 

(Martinet  1992). 

Figure  2.16  shows  the  step  well  structure  and  the  calculated  energy  levels  using 
the  Matlab  program  developed  using  the  “staircase”  approximation  for  the  same  structure 
used  by  Martinet  et  al.  (Martinet  1992). 
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Quantum  Well  without  External  Electric  Field 
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Figure  2. 16  Calculated  energy  levels  of  the  Al^^Ga^^^As  /  Al^^^Ga^^f^As  IGaAs  step 

quantum  well. 


The  calculated  value  of  E^-E^  of  107  meV is  close  to  the  zero- field  absorption 
peak  of  115  meV  observed  experimentally.  The  experimentally  observed  shift  of  the 
absorption  peak  (Martinet,  1992)  for  different  electric  fields  in  the  step  quantum  well  is 
shown  in  Figure  2.17. 


27 


Photoresponse  (normalized) 


Energy  (meV) 

Figure  2.17  Normalized  responsitivity  of  the  asymmetric  step  quantum  well  under 
different  positive  and  negative  biases.  The  peak  shift  is  nearly  four  times  the  spectral 
width.  Zero  field  absorption  (bold  line)  is  included  as  a  reference,  (Martinet ,  1992). 


Normalized  Program  Output  for  Different  Electric  Flelde 


Figure  2.18  Plot  of  bound- to-bound,  1— >2  normalized  absorption  transitions  with 

F  =  20meV  done  with  the  Matlab  program. 
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Figure  2.18  shows  the  calculated  absorption  coefficient  (normalized)  for  the  same 
set  of  bias  voltages  as  in  the  experiment.  The  experiment  shows  that  the  peak  position 
changes  from  93  meV  to  145  meV  when  the  bias  is  changed  from  42  kV  I  cm  to 
-40  kV  I  cm .  The  calculated  value  of  the  peak  changes  from  86  meV  to  125  meV  for 
the  same  bias  values.  The  experimental  and  calculated  peak  shifts  amount  to  52  meV  and 
39  meV,  respectively.  The  difference  between  the  two  values  may  be  due  to  the 

comparison  of  the  photocurrent  versus  the  absorption  coefficient.  Nonetheless,  the  results 
show  that  the  model  predictions  are  in  reasonably  good  agreement  with  the  experimental 
observation,  indicating  the  accuracy  of  the  calculation. 

H.  SUMMARY 

In  this  chapter  the  transfer  matrix  method  was  derived  for  an  arbitrary  potential 
well.  This  technique  was  employed  to  calculate  the  energy  levels  and  wavefunctions  of  a 
step  potential  well.  Intersubband  transition  selection  rules  in  both  symmetric  and 
asymmetric  potential  wells  and  how  they  affect  the  design  of  the  detector  were  also 
discussed.  Another  significant  design  parameter,  the  absorption  coefficient,  was 
calculated  for  both  symmetric  and  asymmetric  wells.  External  bias  and  its  effects  on  the 
quantum  wells,  energy  levels,  and  absorption  coefficients  were  discussed.  The  use  of 
“staircase”  approach  to  calculate  the  absorption  coefficient  is  found  to  be  in  good 
agreement  with  the  experimental  observation.  In  the  following  chapter,  the  design  of  the 
tunable  bandwidth  quantum  well  detector  will  be  discussed. 
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m.  DESIGN  OF  A  TUNABLE  BANDWIDTH  AlGaAs/GaAs  STEP 
QUANTUM  WELL  INFRARED  PHOTODETECTOR 


A.  INTRODUCTION 

As  discussed  in  Chapter  II,  intersubband  absorption  results  from  transitions 
between  the  energy  levels  of  a  quantum  well  formed  within  the  conduction  or  valence 
band.  This  section  presents  detailed  calculations  of  energy  level,  wavefunction  and 
oscillator  strength  in  an  AlGaAsI  GaAs  step  well  structure.  The  use  of  an  electric  field  to 
tune  the  detection  bandwidth  will  also  be  described.  First,  however,  the  calculations  of 
effective  mass  and  barrier  height  as  a  function  of  Al  composition  in  the  Alfiay_^As  layer 
will  be  discussed. 

B.  EFFECTIVE  MASS  AND  BARRIER  HEIGHT  CALCULATION 

Normally  the  band  offsets  of  the  conduction  band  (AEJ  and  valence  band  (AEJ 
of  Alfic\_^As  /  GaAs  quantum  well  structures  can  be  calculated  by  using  the  total  band 
offset  {AE^ ) .  Accepted  values  of  conduction  aid  valence  band  offset  ratios  are  as 
follows  (Meyer,  2001): 

AE^  =  0.6  AE^ 

(3.1) 

AE^  =  0.4AE^. 

The  bandgap  of  Alfia^_^As,E^{x),  can  be  approximately  calculated  using  Vegard’s 
Law  as  (Adachi,  1994) 

(•^)  ~  ^^g(AlAs)  +  (1  ~  ^g(GaA.s)  ’  (3-2) 

where  x  is  the  A/ molar  fraction  and  ^giCaAs)  are  the  bandgaps  of  AZA^  and 

GaAs ,  respectively,  at  the  F  -  point. 

Using  the  experimentally  observed  band  gaps  of  A/As  (2.671  eV)  and 
GaAs  (1.424  eV) ,  the  band  gap  of  Alfia^_^As  at  the  F  -  point  is  given  by  (Singh,  1993) 

F^(x)  =  1.424  +  1. 247 X  eV  (3.3) 
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The  graph  in  Figure  3.1  shows  the  band  gap  of  the  Alfia^_^As  layer  as  a  function  of 
Al  composition  where  the  composition  x  changes  from  zero  to  one  hundred  percent. 


Energy  Gap  vs.  Molar  FracJion  of  fit 


Figure  3.1  Energy  gap  at  F  -  point  as  a  function  of  Al. 

1.  Effective  Mass  Calculation 


The  effective  mass  of  the  Alfia^_^As  alloy  can  also  be  found  using  Vegard’s 

Law  as 


1 

* 

^(AlGoAs  ) 


X 


^(AlAs) 


+ 


1-x 

^(GaAs) 


(3.4) 


The  inverse  relationship  is  due  to  the  fact  that  the  effective  mass  and  energy  are  inversely 
related,  where 


^(AlAs) 


=  0.15m^ 


=0.067m^ 


(3.5) 
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are  the  F  point  effective  masses  of  AlAs  and  GaAs ,  respectively(Singh,  1993). 
Figure  3.2  shows  the  plot  of  effective  mass  as  a  function  of  Al  composition. 


Effective  Mass  vs  Molar  Fraction  of  Al 


Figure  3.2  Effective  mass  vs.  molar  fraction  of  A/. 
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IV.  FINAL  DESIGN  OF  THE  STRUCTURE 


In  the  following,  we  will  diseuss  the  design  of  a  tunable  bandwidth  detector  with 
center  wavelength  near  8.5  [im  .  The  most  convenient  material  system  is  AlGaAs !  GoAs 
due  to  their  lattice  match  and  to  mature  growth  technology.  In  order  to  demonstrate  the 
band  tuning  capability,  a  bound-to-bound  transition  in  a  step  quantum  well  will  be  used. 
Figure  4.1  shows  a  step  well  with  peak  absorption  near  8.9  \im  .  This  structure  was 
optimized  to  achieve  a  relatively  large  absorption  coefficient  with  peak  absorption 
between  8  and  9  \im  .  The  first  excited  state  lies  about  24  meV  below  the  barrier  edge. 

Quantun  Wei  without  External  Electric  Field 
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Figure  4.1  Optimized  quantum  well  showing  the  energy  levels,  centroids,  and 

lambda. 

Figure  4.2  shows  the  calculated  absorption  coefficient  of  the  structure  for  a  set  of  positive 
and  negative  bias  voltages. 
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Bourx)  to  Bound  Absorption  With  Applied  Bias 


Figure  4.2  Absorption  coefficient  for  bound- to-bound  intersubband  transition  for  the 

step  well  shown  in  Figure  4.1  for  a  set  of  bias  voltages. 


Table  4.1  shows  E^,  and  oscillator  strength,  /  ,  for  the  quantum  well  structure  given 
in  Figure  4.1  under  various  applied  biases. 


Energy  states 
and  oscillator 
strength 

Zero 

Biased  Well 

IE  /  |lm 

2V  /  |lm 

3E  /  |lm 

V 

-V 

V 

-V 

V 

-V 

E,{ey) 

0.0766 

0.0833 

0.0767 

0.0926 

0.0792 

0.1017 

0.0818 

E,{eV) 

0.2160 

0.2184 

0.2220 

0.2211 

0.2275 

0.2224 

0.2334 

f 

0.1343 

0.1146 

0.1523 

0.0960 

0.1831 

0.0766 

0.1798 

Table  4.1  Energy  states  and  oscillator  strengths  under  different  biases  for  the 
quantum  well  structure  shown  in  Figure  4.1. 
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It  can  be  seen  in  Figure  4.2  that  the  absorption  peak  shifts  with  applied  bias  due  to 
the  linear  Stark  effeet  in  the  step  quantum  well.  However,  the  bandwidth  of  the 
absorption  peak  remains  nearly  the  same.  In  order  to  aehieve  tunable  bandwidth,  a  two 
step  quantum  well  was  arranged  as  shown  in  Figure  4.3. 


Quartum  Well  v»ithou  Extefnal  Electric  Field 


Figure  4.3  Final  design  of  the  asymmetrie  quantum  well  strueture. 


In  this  case,  under  a  given  bias  the  absorption  peaks  of  the  two  wells  move  in 
opposite  direetions.  Figure  4.4  shows  the  calculated  absorption  coefficient  of  the 
combined  structure  as  a  function  of  photon  energy  for  a  set  of  bias  voltages. 
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Boutxl  to  Bomd  Absorption  With  Applied  Bos 


Figure  4.4  Calculated  absorption  coefficients  for  the  combined  step  quantum  well. 


It  can  clearly  be  seen  that  the  bandwidth  increases  as  the  bias  increases.  At  10 
kVIcm,  the  bandwidth  is  nearly  doubled.  For  higher  electric  fields,  the  absorption 
coefficient  near  the  center  drops  due  to  a  large  shift  of  the  two  peaks.  The  details  of  the 
layer  structure  needed  to  fabricate  the  tunable  bandwidth  quantum  well  detector  are 
shown  in  Table  4.2. 
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Layer 

and  Thickness 

Molar  Fraction 

Layer 

m*  (kg) 

Eg  (eV) 

AlGaAs  (barrier) 
5nm 

0.32 

Alo.32Gao.68As 

0.08 14mo 

1.8230 

AlGaAs  (step) 

7  nm 

0.25 

Alo.25Gao.75As 

0.0778mo 

1.7357 

GaAs 

4.8  nm 

-- 

GaAs 

0.067mo 

1.4240 

AlGaAs  (barrier) 

5  nm 

0.32 

Alo.32Gao.68As 

0.08 14mo 

1.8230 

GaAs 

4.8  nm 

-- 

GaAs 

0.067mo 

1.4240 

AlGaAs  (step) 

7  nm 

0.25 

Alo.25Gao.75As 

0.0778mo 

1.7357 

AlGaAs  (barrier) 
5nm 

0.32 

Alo.32Gao.68As 

0.08 14mo 

1.8230 

Table  4.2  The  layered  structure  of  a  tunable  bandwid 

ith  step  quantum  well  detector. 

In  summary,  the  design  of  a  step  quantum  well  structure  for  tunable  bandwidth 
detection  was  introduced.  The  tunability  of  the  detector  by  applying  a  bias  to  the  quantum 
structure  was  also  demonstrated. 
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V.  CONCLUSION 


The  transfer  matrix  method  was  derived,  and  energy  levels  and  wavefunctions  for 
an  asymmetric  arbitrary  potential  well  were  calculated.  Intersubband  transition  selection 
rules  and  their  effect  on  the  design  process  were  also  discussed.  The  absorption 
coefficient,  which  is  also  important  in  the  design  of  a  quantum  well,  was  calculated  for 
both  symmetric  and  asymmetric  wells.  The  “staircase”  potential  method  was  used  to 
present  both  the  effects  of  an  external  bias  and  the  shift  in  the  detection  wavelength  due 
to  the  various  applied  external  electric  fields  across  the  detector.  External  bias  and  its 
effects  on  the  quantum  wells,  energy  levels,  and  absorption  coefficients  were  explained  in 
detail.  The  validity  of  the  program  was  demonstrated  through  comparison  with 
experimental  data  and  the  best  possible  design  of  the  step  quantum  well  structure  for 
broadband  detection  was  found. 

A  quantum  well  structure  having  a  broadband  detection  capability  was  designed 
using  two  step  quantum  well  structures  arranged  opposite  to  each  other.  This  kind  of 
structure  allowed  both  blue  and  red  Stark  shifts  of  the  transition  energy  producing  a 
tunable  bandwidth.  The  calculated  absorption  spectral  bandwidth  as  a  function  of  the 
applied  bias  shows  that  the  peak  width  can  be  varied  from  1.27  \im  to  2.17  \im  ,  which 

is  twice  the  bandwidth  at  zero  bias. 

Finally,  the  analysis  indicates  the  possibility  of  fabricating  a  tunable  bandwidth 
infrared  detector  using  asymmetric  quantum  wells. 
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APPENDIX.  [SELECTED  PROGRAMS] 


%THIS  PROGRAM  FINDS  BOUND  ENERGY  LEVELS  IN  ANY  KIND  OE  STEP  WELL 
%by  ATAKAN  KONUKBAY  ©  2002 
clear  all,clc,close  all; 
global  N  X  El  h  Bin  v 

N=input('Enter  the  number  layers>  ’);  %ENTER  1  EOR  SQUARE, 2  EOR  STEP  QUANTUM  WELL 

e=1.602e-19;  %Coulomb 
h=197.33e-9;  %Planck's  Constant 

%MOLAR  RATIOS  (ENTER  VALUES  AS  .35) 

XX(I)=input('Enter  the  percentage  of  A1  in  region  I  >'); 

XX(2)=input('Enter  the  percentage  of  A1  in  region  II  >'); 

XX(3)=input('Enter  the  percentage  of  A1  in  region  IV  >'); 

%THE  ENERGY  GAP  IN  EACH  REGION  OE  THE  ASYMMETRIC  QUANTUM  WELL 

Eg  I = 1 .424+ 1 . 247  *XX(  I ) ;  %  AlxGa(  I- x)  As 

Eg2=I. 424+1. 247*XX(2);  %AlxGa(I-x)As 

Eg3=1.424;  %GaAs 

Eg4=  1 .424+ 1 . 247  *XX(  I ) ;  %  AlxGa(  I-  x)  As 

%EEEECTIVE  MASSES 
M=.5IIe6; 

m(  I)=I/((XX(  I)/(0. 15*M))+((  I-XX(  I))/(0.067*M))); 
m(2)=I/((XX(2)/(0.I5*M))+((I-XX(2))/(0.067*M))); 
m(3)=0.067*M; 

m(4)=I/((XX(3)/(0.I5*M))+((I-XX(3))/(0.067*M))); 

%PQTENTIAL  IN  CONDUCTION  BAND 
v(I)=0.6*(EgI-Eg3); 

v(2)=0.6*(Eg2-Eg3);  %Set  origin  @  0 
v(3)=0;  %.6*(Eg3-Eg2); 

v(4)=0.6*(Egl-Eg3); 


a=input('Enter  the  width  from  0->a  in  Angstrom>  '); 
b=input('Enter  the  width  from  a->b  in  Angstrom>  '); 

x=[0  a  a+b]*le-10; 


E=0.0001; 

Ein=input('Enter  the  energy  (eV)>  '); 

Efin=input('Enter  the  external  E  field  (le6  V/m)>  ');  %ENTER  AS  1,2... ETC. 
El_=Ein*le4; 
while  E<=Ein; 
for  d=l:El_ 

T=l; 

for  n  =  1:N+1 

k(n)=sqrt((2*m(n))*(E-v(n))/((h''2))); 
k(n+ 1  )=sqrt((2  *m(n+ 1 ))  *  (E-v(n+ 1  ))/((h''2))) ; 
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A(n)=(exp(i*(x(n)*(k(n)-k(n+l)))))*((k(n+l)*(m(n)))+(k(n))*(m(n+l)))/(2*(k(n+l))*(m(n))); 

B(n)=(exp(-(i*(x(n)*(k(n+l)+k(n))))))*((k(n+l)*(m(n)))-(k(n))*(m(n+l)))/(2*(k(n+l))*(m(n))); 

C(n)=(exp(i*(x(n)*(k(n)+k(n+l)))))*((k(n+l)*(m(n)))-(k(n))*(m(n+l)))/(2*(k(n+l))*(m(n))); 

D(n)=(exp(-(i*(x(n)*(k(n)-k(n+l))))))*((k(n+l)*(m(n)))+(k(n))*(m(n+l)))/(2*(k(n+l))*(m(n))); 


MM(n)={[A(n)  B(n);C(n)  D(n)]}; 
end 


ZZ=1;  %Unit  matrix 
for  n=l:N+l; 

ZZ=MM{n}*ZZ; 

end 

M22=real(ZZ(2,2)); 

TT(d)={[T*M22,E]}; 

Tt(d,:)=TT{d}; 

E=E+.0001; 

end 

end 


zz=l; 

for  z=l:El_-l 
ttl(z)=Tt(z,l); 
tt2(z+l)=Tt(z+l,l); 
ttm=ttl(z)*tt2(z+l); 
if  ttm<0 
E(zz)=[Tt(z,2)]; 
zz=zz+l; 

end 

end 

end 

clc 

MLDRAWEINE 

%MLDRAWEINE 

%THIS  PROGRAM  DRAWS  THE  SHAPE  OE  THE  QUANTUM  WELL  &  DIVIDES  THE  WIDTH  OE 
THE  STRUCTURE  INTO  MANY  STEPS 
%by  ATAKAN  KONUKBAY  ©  2002 

clc; 

e=  1 .602e  -1 9;%Coulomb 
Edp=Efin;  %Bias 
delta=Ie-I0; 

X{I,I}=(x(I)-50e-I0):delta:x(I); 
for  n=2:N+I 

X{  I,n}=x(n-I):delta:x(n); 
end 

X{  I  ,n+ 1 }  =x(  I  ,n):delta:  (x(  I  ,n)+50e  - 10); 
for  n=I:N+2 
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XL(l,n)=length(X{l,n}); 

Y{  l,n}=(Edp*X{  l,n})*le6;%cnst; 
YL(l,n)=length(Y{l,n}); 
end 

for  n=l:N+2 

lgth=length(  Y  { n } ) ; 
Y_{l,n}=[Y{l,n}(l:lgth)]+v(n); 

end 


x_=[]; 
y_=[]; 
for  n=l:N+2 

x_=union(x_,X  { n } ); 
y _ =cat(2,y _ ,  Y_  { n } ) ; 

end 

for  n=l:N+2 
hold  on; 

plot(X{  l,n},(Y{  l,n}+v(l,n))); 
end 


for  n=2:N+2 
LXf=XL(l,n); 

LYf=YL(l,n); 
hold  on 

line([X{n-l}(l,XL(l,n-l)),X{n-l}(l,XL(l,n-l))],[Y{n-l}(l,YL(l,n-l))+v(n-l),Y{n}(l,l)+v(n)]) 
hold  on 
end 

for  n=l:length(E) 
switch  N 
case  1 
hold  on 

line([x(l),x(2)],[E(n),E(n)],'ColorVr') 
hold  on 
case  2 
ifE(n)<v(2) 
hold  on 

line([x(2),x(3)],[E(n),E(n)], 'Color', V) 
hold  on 
else 

line([x(l),x(3)],[E(n),E(n)], 'Color', 'r') 

end 

end 

end 

title(sprintf('Quantum  Well  With  External  Electric  Eield  %g  10^6  V/m',Efin)); 

for  n=l:N+2 
QA{n}=Y{n}+v(n); 

L(n)=length(Q  A  { n } ) ; 
end 
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switch  N 


case  1 

AQ=QA{2}; 

AQ(1)=[]; 

AQ((L(2>  !))=[]; 

QA{2}=AQ; 
v_=[]; 
case  2 

AQ=QA{2}; 

AQ(1)=[]; 

QA{2}=AQ; 

AQ=QA{3}; 

AQ(1)=[]; 

QA{3}=AQ; 

v_=[]; 

end 

for  n=l:N+2; 

v_=cat(2,v_,QA{n}); 

end 

L=length(v_); 

vl=v(l); 

xlabel(sprintf('a  =  %g  nm  ,  b  =  %g  nm',x(2)*le9,(x(3)-x(2))*le9)); 
ylabel(sprintf('Energy  (eV)',Efin)); 
set(findobj(gca, 'Type', 'line', 'Color', [0  0  1]),... 

'Color','blue',... 

'LineWidth',1) 

text(x(3)+28e-10,0+.01,sprintf('%g  /  %g  /  %g',XX(l),XX(2),XX(3)),... 

'Horizontal  Alignment' ,  'left',  'EontS  ize' ,  8) 
for  n=l:length(E) 

text(x(3)+le-10,E(n),sprintf('E%g=%g  ev',n,E(n)),... 

'Horizontal  Alignment' ,  'left',  'EontS  ize' ,  8) 
end 

grid  off; 

%THIS  PROGRAM  CALCULATES  THE  ENERGY  LEVELS  EOR  MULTI  STEP  WELL 
%by  ATAKAN  KONUKBAY  ©  2002 

global  N  V  X  El  h  Lin 
N=length(x_)-1; 


e=1.602e-19;  %Coulomb 
h=197.33e-9;  %Planck's  constant 

x=x_; 

%******mass 

m_l(l:XL(l))=m(l); 
m_2(l:(XL(2)-l))=m(2); 
m_3(l:(XL(3)-l))=m(3); 
m_4(l:(XL(4)))=m(4); 
mm=[m_l  m_2  m_3  m_4]; 
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v=v_; 


E=0.0001; 

Edrw=0; 

Ein=input('Enter  the  energy  level  (eV)>  '); 

Efin=0 
El_=Ein*le4; 
while  E<=Ein; 
for  d=l:El_ 

T=l; 

for  n  =  1:N+1 

k(n)=sqrt((2*mm(n))*(E-v(n))/((h''2))); 
k(n+ 1  )=sqrt((2  *mm(n+ 1 ))  *  (E-v(n+ 1  ))/((h''2))) ; 

A(n)=(exp(i*(x(n)*(k(n)-k(n+l)))))*((k(n+l)*(mm(n)))+(k(n))*(mm(n+l)))/(2*(k(n+l))*(mm(n))); 

B(n)=(exp(-(i*(x(n)*(k(n+l)+k(n))))))*((k(n+l)*(mm(n)))- 

(k(n))*(mm(n+l)))/(2*(k(n+l))*(mm(n))); 

C(n)=(exp(i*(x(n)*(k(n)+k(n+l)))))*((k(n+l)*(mm(n)))-(k(n))*(mm(n+l)))/(2*(k(n+l))*(mm(n))); 

D(n)=(exp(-(i*(x(n)*(k(n)- 

k(n+l))))))*((k(n+l)*(mm(n)))+(k(n))*(mm(n+l)))/(2*(k(n+l))*(mm(n))); 

MM(n)={[A(n)  B(n);C(n)  D(n)]}; 
end 

ZZ=1;  %Unit  matrix 
for  n=l:N+l; 

ZZ=MM{n}*ZZ; 

end 

M22=real(ZZ(2,2)); 

TT(d)={[T*M22,E]}; 

Tt(d,:)=TT{d}; 

E=E+.0001; 

end 

end 


zz=l; 


for  z=l:El_-l 
ttl(z)=Tt(z,l); 
tt2(z+l)=Tt(z+l,l); 
ttm=ttl(z)*tt2(z+l); 
if  ttm<0 
E(zz)=[Tt(z,2)]; 
Edrw(zz)=E(zz)+.00 1 ; 

zz=zz+l; 

end 

end 
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%THIS  PROGRAM  CALCULATES  THE  ALPHA  EOR  BOUND  TO  CONTINUUM  TRANSITIONS 
%by  ATAKAN  KONUKBAY  ©  2002 


clc, close 

e=L602e-I9;  %C 

Nd=IeI8/Ie-6;  %I/m3 

hbar=L055e-34;  %J.s 

EoE=8.85e-I2*sqrt(I3);  %P/m 
c=3e8  %m/s 


ENG=v(I)+.000I:.000I:v(I)+.02;%-(E(2)-E(I)):.00I:2*((E(2)-E(I))); 

Alpha=Nd*e''2*hbar''2*sqrt((2*m(I)*9.IIe-3I/.5IIe6)./((((ENG-E(I)).*e).''2).*((ENG- 

v(I)).*e)))./(2*((m(3)*9.IIe-3I/.5IIe6)''2)*EoE*c)*AvXRegBCA2 

plot(ENG,Alpha/Ie2) 
grid  on 

xlabel('h\prime\omega  [eV]') 
ylabel('Absorption  [cm''{  -I }]') 
title('Bound  to  Continuum  Absorption') 

%THIS  PROGRAM  CALCULATES  THE  ALPHA  EOR  BOUND  TO  BOUND  TRANSITIONS 
%by  ATAKAN  KONUKB  A  Y  ©  2002 

clc, close 

Nd=lel8/le-6  %l/m3 
e=1.602e-19;  %C 
hbar=1.055e-34  %J.s 
T=20e-3*e/2  %J 
EoE=8.85e-12*sqrt(13);  %E/m 
c=3e8  %m/s 

HW=E(  l)+.000 1 :  .000 1  :E(2); 

G=T./((((E(2>E(l)-HW)*e).''2)+((T)''2)); 

CB_B=(Nd*(e''2)*(hbar^3))/(((m(3)*9.11e-31/.511e6)''2)*EoE*c); 

Alpha=(CB_B.*G.*AvXRegMB_B''2)./(HW.*e) 

plot(HW,Alpha/le2) 
grid  on 

xlabel('h\prime\omega  [eV]') 
ylabel('Absorption  [cm''{  -1 }]') 
title('Bound  to  Bound  Absorption') 

%THIS  PROGRAM  CALCULATES  THE  WAVEEUNCTION  EOR  MULTISTEP  WELL 
%by  ATAKAN  KONUKBAY  ©  2002 

clc; 

figure; 

delta=le-10; 


trans=MM; 
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trs=MM; 
for  n=2:N+l 

trans{  l,n}=trs{  l,n}*trans{  l,n-l}; 
end 

A(1,1)=0; 

B(l,l)=l; 
for  n=l:N+l 
wavtrans=trans  { l,n}; 

A(  1  ,n+ 1  )=wavtrans(  1 ,2)  *B  ( 1 , 1 ) ; 

B(  l,n+  l)=wavtrans(2,2)*B(  1 ,1); 
end 

A(1,1)=0; 

B(l,l)=l; 

B(  1  ,(XL(  1  )+XL(2)+XL(3))  - 1  :n+ 1  )=0; 

X{l,l}=-x(2):delta:x(l); 
for  n=2:N+l 

X{  l,n}=x(n-l):delta:x(n); 
end 

X{  l,n+l  }=x(l,n):delta:(x(l,n)+(x(l,n)-x(l,n-l))); 

TotalIntegral=0; 
for  n=l:N+2 

W{  l,n}=A(l,n)*exp(i*k(l,n)*X{  l,n})+B(l,n)*exp(-(i*k(l,n)*X{  l,n})); 

%%%***for  the  wave  used  in  calculating  the  oscillation  with-B  parameter 
W_BC{  l,n}=A(l,n)*exp(i*k(l,n)*X{  l,n})-B(l,n)*exp(-(i*k(l,n)*X{  l,n})); 

%%%*** 


W_{  l,n}=conj(W{  l,n}); 

integr{  l,n}=(W{  l,n}.*W_{  l,n})*delta; 

integregn{  l,n}=sum(integr{  l,n}); 

T  otalIntegral=integregn  { 1  ,n } +T  otalintegral ; 
end 

bl  =  l/(sqrt(TotalIntegral)); 

normwave=0; 

for  n=l:N+2; 

WavReal{  l,n}=bl*W{  l,n}; 

%%%***Now  multiplying  with  normalization  factor  to  get  the  wave 
WavRealBC  { 1  ,n }  =b  1  *W_BC  { 1  ,n } ; 

%%%*** 


WR_{  l,n}=conj(WavReal{  l,n}); 

Integ{  l,n}=(WavReal{  l,n}.*WR_{  l,n})*delta; 
Integral!  l,n}=sum(Integ{  l,n}); 
normwave=Integral  { 1  ,n  }+normwave; 
end 

for  n=l:N+2 
hold  on; 

plot(X{  l,n},WavReal{  l,n}); 

hold  on; 

end 
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